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Abstract 

The mixed phase of layered superconductors with no magnetic screening is studied through a 
partial duality analysis of the corresponding XY model in the presence of random pinning centers. 
Sufncently weak Josephson coupling between adjacent layers results in an entangled vortex solid 
that exhibits weak superconductivity across layers. The corresponding vortex liquid state shows 
an inverted specific heat anomaly that is a precursor to the weak superconductor. 
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INTRODUCTION 



The possibility that the vortex lattice of a pristine layered superconductor with external 
magnetic field oriented perpendicular to the layers decouples into an incoherent stack of 
two-dimensional (2D) vortex lattices before it melts has been proposed many times in the 
literature Monte Carlo simulations and a duality analysis of the corresponding 

layered XY model with uniform frustration that describes the extreme type-II limit find 
that the decoupled vortex lattice state does not exist in practice however. The presence 
of random pinning centers [fj ?1 result in quenched-in entanglements of the vortex lines that 
can reverse this conclusion [8]. 

In this paper, we shall show that random point pins drive a dimensional cross-over tran- 
sition in the zero-temperature limit between defective vortex lattice states that exhibit weak 
versus strong superconductivity across layers. The magnetic field at the dimensional cross- 
over transition notably decreases with an increase in the level of random point pinning. This 
is consistent with recent experimental determinations of how the "fishtail" feature that is 
shown by the irreversible magnetization of organic superconductors in magnetic field ori- 
ented perpendicular to the layers moves with a change in the level of material disorder jjj. 
We also find that the introduction of correlated pins can act to increase the dimensional 
cross-over field in the zero-temperature limit. A similar effect has been recently shown by 
the dimensional cross-over field extracted from flux creep dynamics in high-temperature 

rn 

superconductors [10|. Last, we find that the zero-temperature cross-over in magnetic field is 
inextricably linked to the existence of an inverted specific heat anomaly within the vortex 

n 

liquid phase j8|. The latter accounts for the experimental observation of a similar thermal 



anomaly inside the vortex liquid phase of the high-temperature superconductor YBCO 
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THEORY 



We shall describe the mixed phase of layered superconductors with random pins by the 
corresponding XY model with uniform frustration over the cubic lattice. Both fluctuations 
of the magnetic induction and of the magnitude of the superconducting order parameter are 
neglected within this approximation. The model is therefore valid deep inside the interior 
of the mixed phase. The thermodynamics of the three-dimensional (3D) XY model with 
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anisotropy and with uniform frustration is determined by the superfluid kinetic energy 

Efy = - E E-WA^-^]| n (1) 

fi=x,y,z r 

which is a functional of the superconducting phase 0(r) over the cubic lattice. The local 
phase rigidities within layers, J x and J y , are assumed to be constant over most of the 
nearest-neighbor links, with the exception of those links in the vicinity of a pinning site, 
while the Josephson coupling J z on links between adjacent layers is assumed to be constant 
and weak. The vector potential = (0, 2irfx/a, — b»x) represents the magnetic induction, 
with a component B± = &of /a 2 oriented perpendicular to the layers, and with a component 
B\\ = ($ /2vrc/)6|| oriented parallel to the layers. Here a denotes the square lattice constant 
for each layer, which is of order the zero-temperature coherence length, while d represents 
the separation between layers. Also, $o denotes the flux quantum, and / denotes the 
concentration of planar vortices per site. 

Partial Duality Analysis. The layered XY model (0) can be effectively analyzed in the 
selective high-temperature limit, ksT 3> J z , through a partial duality transformation j^j]. 
This leads to a dilute Coulomb gas (CG) ensemble that describes the nature of the Josephson 
coupling in terms of dual charges that live on the vertical links. In particular, thermodynamic 
quantities are obtained by taking appropriate derivatives of the logarithm of the partition 
function 

zp y [o] = [h{J z /k B T)f ■ Z CG [0] ■ n,z&[o] (2) 

for the layered superconductor, expressed here as a product of the partition function for a 
layered CG ensemble 0] 

Zcg[0]= E ^•nAfel-e-E^A, (3) 

{n.(r)} 

(2) 

with the contributions from each layer I in isolation, Z X y, along with a modified Bessel 
function, Iq{x), raised to the total number of links between adjacent layers, M' . Above, the 
dual charge, n z (r, I), is an integer field that lives on links between adjacent layers / and / + 1 
located at 2D points r. The CG ensemble is weighted by a product of phase auto-correlation 
functions 

Q[qi] = (exp[i^g z (r)0(f,O]} J=o (4) 
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for isolated layers / probed at the dual charge that accumulates onto that layer: 



qi{r) = n z (r,l - 1) -n z (r,l). (5) 

It is also weighted by a bare fugacity yo that is raised to the power iV[n z ] equal to the total 
number of dual charges, n z = ±1. The fugacity is given by yo = J z /2ksT in the selective 
high-temperature regime, J z <C &bT, reached at large model anisotropy. 

Defective Vortex Lattices. A recent analysis of the 2D limit of the XY model (jT|) . J v = 0. 
finds that a phase coherent vortex lattice state with quenched- in unbound dislocations jl^ 
exists in the zero-temperature limit if the random pinning is weak enough 8]. This defective 
vortex lattice state exhibits 2D phase auto-correlation functions (jlj) of the form Ci[q] = 
\Ci[q]\exp[iJ2fq(r)(po(r, I)], where the magnitude decays algebraicly like 

\Ci[q} \ =9o + -exp[r]2D ffMMW^) ?( 2 )] ( 6 ) 

(1.2) 

asymptotically, and where 0o(r, Z) represents the zero-temperature configuration of layer / 
in isolation. Above, r] 2 D is the correlation exponent that vanishes linearly with temperature 
in the zero-temperature limit, and r$ is the natural ultraviolet scale. Also , go is equal to the 
ratio of the 2D phase rigidity with its value at zero temperature, while n + is equal to half 
the number of probes in q(r). At this stage it becomes useful to re-express the layered CG 
enesemble (JHJ) by replacing Ci[q] with its magnitude (jBJ), and by compensating this change 
with the additional replacement of —A z (f,l) with 

= I + 1) - Mr, I) ~ A z (f, I). (7) 

A Hubbard-Stratonovich transformation of the CG partition function (J3J) then reveals 
that it is equal to the partition function for a renormalized Lawrence-Doniach (LD) model, 
Zld[0] = / V9 e~ ELD ^ kBT , up to a factor that is independent of the Josephson coupling, J z . 
The LD energy functional is given by(j| 



(8) 



where 9i t i+i = <pfj +1 + 9i + i — 9i. The 2D phase rigididity above is related to the 2D correlation 
exponent by p^ D > = ksT /2nr]2D- Also, A = 7' a is the Josephson penetration depth, where 
7' = {J/J z ) l l 2 is the model anisotropy parameter with respect to the macroscopic phase 
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rigidity of an isolated layer at zero-temperature, J. A standard thermodynamic analysis of 
the product of partition functions O) then yields that the strength of the local Josephson 



coupling is given by (cos0 M+ i) = y + # o (cos0 M+ i). Here 

fo+i(r) = 4>(f, l + l)- <!>(?, I) - Mr, I) (9) 

is the gauge- invariant phase difference across adjacent layers, while the overbar notation 
represents a bulk average. It can also be shown|3] that the phase rigidity across a macroscopic 



number of layers is equal to pjr/J z = g ( cos $M+i)- 



To compute the bulk average (cos 6^+1) at low temperature, we must first determine the 
configuration df'\r) that optimizes -Eld- Eq. (JHJ) implies that it satisfies the field equation 

- V^ 0) + A - 2 [sin ef\, - sin J = 0. (10) 

In the weak-coupling limit, A — > 00, we therefore have that 6\ \ (r) is constant inside of 
each layer. Also, by analogy with 2D melting phvsics the presence of quenched-in 

unbound dislocations in the pinned vortex lattice [8j]|l2] implies that the quenched-in inter- 
layer order parameter expi<fiii +1 shows only short-range correlations on average over the 
bulk. Specifically, we have 



exp^g^l)] ■ exp[-#; ) ( ; ) +1 (2)] = e »i«" e~ r ^ (11) 

asymptotically, where 1$ is a 2D disorder scale set by the density of unbound dislocations 
between adjacent layers (see next section). The absence of long-range order along adjacent 
layers on average over the bulk then implies that cos 9j°i + x — {56^ — S@i+i) • sin 0j at zero 
temperature in the weak-coupling limit [6], where 80i is the deviation with respect to the 
constant background. This yields the perturbative result 

cos<ti=Ao 2 / d 2 r 12 ex P [i0{? +1 (l)] • exp[-i^ 1 (2)] ■ G^(r 12 ) (12) 

at weak coupling by Eq. (|10|) . where = — V~ 2 is the 2D Greens function. After taking 
the separation at which vanishes to be of order Ao, we then obtain the final perturbative 
result 3] 

costf^ - [A - A 2 (6||^) 2 ](/ /A o ) 2 ln(A o // ) (13) 

for the bulk average at zero temperature in the weak-coupling limit, and at weak parallel 
field B\\ <C $o/27r/^(i. Here A and A 2 are numerical constants of order unity: e.g., A = 1 



and A 2 = 3/2 if Eq. (fTTj) is generally valid. By the discussion following Eq. (JUJ), we 
conclude that random point pins result in a vortex glass that exhibits weak superconduc- 

Bfl 
, Pg <C J z , at sufficiently high layer anisotropy 8], A ^> 1$. Notice, 

however, that the quadratic dependence with small parallel field shown by the interlayer 

LD cosine (|13|) implies the absence of line tension for Josephson vortices! We believe that 

this reflects the vortex-glass nature of the phase coherence 2]: (e l< ^. ! +i(°) e~ l ^ l - l+1 ^) = and 

|( e ^u+i(o)e-^u+i(°°))|2 > o. Last, the bulk average at zero temperature cos must pass 

between zero and unity at A ~ l<f> as a function of the anisotropy parameter. This condition 

then defines a decoupling cross-over field 

BUO) ~ (/ /a vx ) 2 ($ o /Ag), (14) 

at which point the reversible magnetization shows a broad diamagnetic peak. Above, a vx = 
a j f 1 / 2 is the planar inter-vortex scale. 

Decoupled Vortex Liquid. Suppose now that the macroscopic phase coherence (jBJ) shown 
by each layer in isolation at low temperature is lost at temperatures above a transition 
temperature T^ 2D \ and that only short-range phase correlations over a scale ^2D exist. A 
high-temperature expansion 1(| of the CG ensemble ® in powers of the fugacity y is then 
possible at temperatures above the 2D glass transition in the weak-coupling limit. At zero 
parallel field in particular, Eqs. (j2J) and (|HJ) yield the expansion 



(cos0 M+1 )= 2/o e^(o,i)-q + i(o,i; 



^a(0,l)-Q +1 (0,l)-Q(2,3)-Q +1 (2,3)- £ Q(0, 1, 2, 3) ■ Q +1 (0, 1, 2, 3) 

1,2,3 1,2,3 

(15) 



for the inter-layer "cosine" to next-leading order in the fugacity. The prime over the sum 
in the last term above excludes double counting of probes (0,1,2,3). Also, next-leading- 
order contributions that couple together three adjacent layers are negligible in the critical 
regime, £ 2 d 3> l<t>, and these are not listed above. Substitution of the form Cj(l,2) = 
g e~ Vl2 ^ 2D e 1 ^ ^ e~ t( ^ ^ for the phase autocorrelation function (J3J into the leading term above 
yields a non-divergent result 



(cos0 M+1 ) ~ gl{J/k B T)[{l^ + ^ 1 )- 1 /A ] 2 , (16) 
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where = £20/2. This implies an anomalous inter-layer contribution to the specific heat 
per volume equal to 

5c$ - 2[1 + {^/l^idln^/dT^j, (17) 

where ej = (cos^+i) ■ J/A^d is the Josephson energy density, and where d denotes the 
spacing in between adjacent layers. It notably shows an inverted specific heat jump that 
is followed by a tailoff at a temperature T p such that £,</>(T p ) ~ if £^ diverges faster than 
[T-TW]- 1 . 

The next-leading-order term for the inter-layer "cosine" ()15j) is negative, divergent, and of 
order — ^g^^/a 6 in the critical region, £^ ^> Z^. This means that it becomes comparable 
to the leading order term (fTBj) at magnetic fields below a cross-over scale 

Bi ~ g {J/k B T)(!&/c&{* /A*) (18) 

that separates 2D from 3D vortex-liquid behavioral The phase correlation length across 
layers grows larger than the separation between adjacent layers at temperatures below this 
point, T x . It then diverges (or jumps to infinity) at a 3D ordering temperature, T g , that lies 
inside of the window [T^ 2D \T X ]. The phase diagram displayed by Fig. 1 summarizes all of 
the above. 



DISCUSSION AND CONCLUSIONS 

The present theory predicts a disorder-driven dimensional crossover between defec- 
tive vortex lattices that show weak versus strong superconductivity across layers. Compari- 
son of the bulk average (fTTJl of the quenched-in phase autocorrelation with the corresponding 
thermal average in the hexatic phase of a pristine vortex liquid [1J] indicates that the 2D 
disorder scale, 1$, that controls the crossover transition is set by the density of unbound 
dislocations between adjacent layers in the vortex lattice. For example, the addition or the 
removal of a straight line of dislocations that traverses all of the layers will not change this 
scale. A useful parameterization of is then 

= W(2a#L), (19) 

where I2D is the correlation length of the quenched-in order parameter exp(i0 o ) 011 average 
over a given layer in isolation, and where x poin t is the concentration of unbound dislocations 
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between adjacent layers relative to the total number of unbound dislocations in adjacent 
layers: e.g., Xpoint = if all of the unbound dislocations form smooth lines across layers, 
while Xp i n t = 1 if sections of such lines are absent. Eqs. (fTlJ) and (fTU|) then yield two 
important predictions: the decoupling crossover field -B^(O) is (i) diminished by an increase 
in the density of unbound dislocations between adjacent layers, l$ 2 , and it (ii) increases 
with an increase in the concentration of unbound dislocations that line-up across adjacent 
layers, 1 — anoint- Prediction (i) is consistent with a recent proposal that the second mag- 
netization peak (or "fishtail") observed in the mixed phase of organic superconductors is 
due to dimensional crossover of the vortex lattice [jj. The second-peak shown by the irre- 
versible magnetization that is measured in this work moves to lower perpendicular magnetic 
field as the concentration of random point pins increases. Also, prediction (ii) is consistent 
with recent experiments that observe a 2D-3D crossover in the flux creep dynamics of high- 
temperature superconductors in magnetic field oriented perpendicular to the layers |lo|. The 
dimensional cross-over field that is extracted experimentally in ref. Q| is found to increase 
after columnar pins are introduced. It is quite possible that these help to line-up a fraction 
of the existing unbound dislocations across layers, which in turn leads to an increase in the 
decoupling field (jl4j) as a result of the increase in the 2D disorder scale 1$. 

The second important result obtained here is the inverted specific-heat anomaly, Eq. (fTTj). 
that is predicted to exist inside of the decoupled vortex liquid phase. Recent experiments on 
2D arrays of Josephson junctions with bond disorder in external magnetic field are consistent 
with a continuous vortex glass transition that has a 2D phase correlation length that diverges 
as ^2D oc [T — Tjj 2D ^]~ U2D at the transition, with a critical exponent i>2D — 2 |l7|. This implies 
an inverted specific heat anomaly, Eq. 1)1 7]). that peaks at £,2d(T p ) = 21$, which is shown 
in Fig. 2. A similar anomaly was observed within the vortex-liquid phase of YBCo[ll|. 
Comparison of the magnitude of that anomaly with Eq. (jl7j) yields a 10% change in the inter- 
layer cosine (fT6|) across the vortex-liquid/vortex-liquid transition j 8]. The above indicates 
that the specific-heat anomaly observed experimentally inside of the vortex-liquid phase of 
YBCO represents not a true phase transition;^, but a crossover instead. 

Noticably absent from the above calculations is the effect of magnetic coupling between 
planar vortices across layers. This in general reduces the entanglement of vortex lines across 
layers, and it will therefore increase the inter-layer phase coherence, Eqs. (I13J1 and (|16|). 
Ref. jp 7 ] shows, however, that the effect of magnetic coupling can be accounted for by 
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a renormalization up in the 2D disorder scale, 1$. Last, a sharp and possibly first-order 
disorder-driven decoupling transition has been found numerically at zero temperature for 
the vortex lattice with magnetic interactions, but with no Josephson coupling^]]. This 
indicates that the disorder- driven dimensional crossover found here in the zero-temperature 
limit survives the addition of weak magnetic coupling. 

We conclude that random point pinning of the vortex lattice in layered superconductors 
results in a 2D-3D crossover transition at low temperature that separates strong from weak 
superconductivity across layers (see Fig. 1). The disorder-driven phenomenon suggests the 
following physical picture: although the tension for a line of dislocations in the defective 
vortex lattice is null on both sides of the crossover, the scale of the roughness in the line is 
large compared to the inter-layer spacing in the 3D phase, while it is of order or less than the 
inter-layer spacing in the quasi-2D phase. Scaling suggests that such quenched-in disorder is 
established first at separations between layers of order l^/j, where 7 = A /d is the physical 
anisotropy parameter. 

The author thanks C. Reichhardt and C. Olson for informative discussions. 
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